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Abstract 

In this work, we study the radiative decay of heavy quarkonium states by using the effective Lagrangian approach. 
Firstly, we construct the spin-breaking terms in the effective Lagrangian for the nP <-> mS transitions and determine 
the some of the coupling constants by fitting the experimental data. Our results show that in Xcj, i/'(25' ), T(25 ), and 
Y(3S ) radiative decays, the spin-breaking effect is so small that can be ignored. Secondly, we investigate the radiative 
decay widths of the cc( \D) states and find the if i^(3770) is a pure 3 Di state its radiative decay into Xcj + J roughly 
preserve the heavy-quark spin symmetry, while if it is a S -D mixing state with mixing angel 12° the heavy quark-spin 
symmetry in its radiative decay and in the radiative decay of i/'(3686) will be largely violated. In the end, we show that 
combining the radiative decay and the light hadron decay of f-wave Xbj(i, 2f) can provide another way to extract the 
O-j information of the color-octet matrix element in the context of non-relativistic QCD (NRQCD) effective theory, and 
q_j our result is consistent with potential NRQCD hypothesis. 
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1. Introduction 

The heavy quarkonium states that are constituted by heavy-quark (Q) and anti-quark (Q) pair provide an ideal 
laboratory to study the dynamics of strong interaction from both perturbative and non-perturbative aspects. In recent 
years, thanks to the large amount data accumulated in electron-positron colliders and hadronic colliders, more accurate 
and new properties of them have been obtained. Especially, many new resonances were discovered, which gave rise to 
a great renewed theoretical interest in studying their spectra and decays (for recent reviews see Ref.[ 1] and references 
therein). 

For the states below open flavor threshold (DD for charmonium and BB for bottomonium), they have relative 
narrow width because they can not decay through the Okubo-Zweig-Iizuka allowed decay mechanism. Their radiative 
decay width could reach hundred kev level, therefore, contributes a considerable branching ratio. On the experimental 
side, the radiative transitions among heavy quarkonia also play an important role in searching for the new states. 
Theoretically, the heavy quarkonium states are approximate nonrelativistic systems. Their annihilation decays are 
sensitive to the wave functions of QQ at small distance, on the contrary the radiative decay can help us to probe the 
behavior of wave functions at long-distance. Besides the intrinsic scale Aqcd, heavy quarkonia are characterized by a 
hierarchy of three energy scales, rtiQ the heavy quark mass, mgvg and fWgVg the typical momentum and energy of the 
heavy quark, where vq <k 1 is velocity of the heavy quark in the rest frame of the heavy meson. The nonrelativistic 
effective field theories, nonrelativistic QCD (NRQCD) (H S 0], and potential NRQCD (pNRQCD) USB] are 
suitable tools to separate the physics in different energy scales. Recently, the magnetic dipole (Ml) transition as 
well as the radiative decay of X(3872) was studied within the framework of pNRQCD[8]. Besides from the model- 
independent perspective, the radiative transitions among heavy quarkonia have already been extensively studied within 
potential model approach (here we refer Ref.|9|] as a comprehensive review). 



Email address: hzgzlhagmail . com (Zhi-Guo He) 



Preprint submitted to elsevier: UB-ECM-PF-70/12, ICCUB-12-097 



February 23, 2012 



In this letter, we will not only calculate the radiative decay widths, but will also do some further analysis by taking 
into account the spin-breaking effect or the S - D mixing effect in the radiative decays. We will also combine the 
radiative decay with the light hadron (LH) decay of the YbjjnP) states to extract the information of the color-octet 
(CO) matrix elements in NRQCD[4] and pNROCD lflol [Till . In this work, we plan to employ the effective Lagrangian 
approach, which can exploit the heavy quark spin symmetry order by order through the expansion of 1/mg. The 
rest of this letter is organized as follows. A brief description of the effective Lagrangian approach will be given in 
Section 2, and then is used to study the S <-» P transitions in cc and bb systems by taking into account spin-breaking 
effect. In section 3, we will study the tfr( } Dj) — > Xcj' + 7 transitions. In section 4, we will relate the transitions of 
Xbj(nP) — > T(nS) + y to the LH decay of Xbj(nP) states to determine the ratios of the CO matrix elements mfHs,{nP) 
to the corresponding color-singlet (CS) matrix elements < H\(nP), where 1-i\(nP) = (xbj(nP)\Oi( 3 Pj)\xbj(nP)), and 
1ri^(nP) - (xi,j(nP)\0%( 3 S \)\xbj(nP)) J4[]. A short summary and conclusion will be presented in the last section. 



2. Effective Lagrangian For Radiative Transitions 

The heavy charmonium states can be classified according to the spectroscopic notation n 2s+l Lj , where n = 1,2,... 
is the radial quantum number, S = 0, 1 is the total spin of the heavy quark pair, L = 0, 1, 2 . . . (or S, P, D . . .) is 
the orbital angular momentum, and J is the total angular momentum. They have parity P = (-1) L+1 and charge 
conjugation C = {-1) L+S . As mentioned in the introduction, NRQCD and pNRQCD are a good starting point to 
describe this system. The LO NRQCD Lagrangian is invariant under S=S U(2)q ® SU(2)q spin symmetry group, an 
approximate symmetry of the heavy quarkonium states, that is inherited in the subsequent effective theories. Hence, 
it is most convenient to introduce hadronic spin-symmetry multiplets, in an analogous way as it was initially done in 



Heavy Quark Effective Theory (HQET) 1 12]. 



For heavy quarkonium states, this formalism was developed in Ref. 111311 . The states have the same radial number 



n and the same orbital momentum L can also be expressed by means of a single multiplet: J^-^ L [ 13], 



(=1 

LV(2L-l)(2L+l)£r 



-1 

1=1 



+ j^-V)^ (i) 

where v M is the four- velocity associated to the multiplet Jfi-f 1 '- (not to be mistaken by vg, the typical velocity 
of the heavy quark in the heavy quarkonium rest frame), K l ^'" flL represents the spin-singlet effective field, and 
ff^'j'^" ,H l ^'" >lL and represent the three spin-triplet effective fields with J - L - 1,L, and L + 1 respec- 

tively. The four tensors are all completely symmetric and traceless and satisfy the transverse condition 

v Ui K£ "« = , r /( //'; '" - 0. (2) 

i = 1, . . . ,L, j = I,..., J. The properties of H and K under parity, charge conjugation and heavy quark spin transfor- 
mations can be easily obtained by assuming that the corresponding transformation rules of the multiplet J^' ■ ■ pL follow 

as: 

(3a) 

j^-kl _^ (-l) i+1 C[7 Ml ... M J r C, (3b) 
±>SJ m .. PL S'\ (3c) 



where C is the charge conjugation matrix (C = iy 2 y° in the Dirac representation), and S e S U(2)q and S ' € S U(2)q 
correspond to the heavy quark and heavy antiquark spin symmetry groups ([S, f] = [S', f] = 0). 

Since we are going to consider the S , P, and D wave states radiative decay,it will be helpful to give the explicit 
expressions of the S-, P-, and D-wave multiplets that follow from Eq.(l). For the L = S case, we have 

J^ l -±l {H ^-K^) l -^L, (4) 



for the L — P case, 



and for L-D case, 



J" = i±t[tify a + -L e ^r Va7l3 H ly + -Ltf - v»)H + K^}^, (5) 



+ [( ^ - v " )H i + w ~ yV)H " ~ lw ~ ^r« H i ] + K 2 V y 5 } ¥ ( 6 ) 

At the leading order of 1 /mg expansion, the radiative transitions between mS and nP states, and between mP and 
nD states can be described by the Lagrangian given in Ref. 1131, 1 1411 : 

£ SP = 2 d"£' mS Tr[J(mS )J l fytP)^ v F n ' + h.c. , (7a) 

m,n 

Z FD = J] 6"°' mP Tr[J a (mP)J afI (nD)]v v F» v + h.c. , (7b) 

where 6"q (Q = c,b) and s mP,nD (Q = c,b) are the coupling constants, and F^ v is the electromagnetic tensor. The 
Lagrangian in Eq. © preserves parity, charge conjugation, gauge invariance and heavy quark and antiquark spin 
symmetry. 

The radiative decays of the low lying S - and P- wave states have been well measured. It will be interesting to do 
some delicate analysis beyond leading order. One important higher order contribution comes from the spin-breaking 
effect which is due to the spin-spin Sq ■ Sg, spin- orbit L • S, and tensor (Se ' r ^ Si} ' r) - i nteract i ons j n pNRQCD 
1 1511 (or in potential models), where Sq and Sg are the spin of the quark and anti-quark respectively, S = Sg + Sq, 
and L is the orbital angular momentum of the heavy meson. To figure out all the spin-breaking terms in the effective 
Lagrangian, it will be more perspicuous to construct them in the rest frame of the heavy meson, where the 4x4 
dimensional space is reduced to the 2x2 dimensional space. In the two component notation, the field J and J 11 is 
simplified as: 

V2 V3 

where & is the Pauli matrix. The Lagrangian in Eq.dTat becomes: 



/ = H ■ & + K , f = (H'± + -^e' 7 X + —H )cr j + K\ , (8) 



£ SP = J] #£" s Ti[J*(mS)J t (nFflE i + h.c. . (9) 

nun 

In the 2x2 dimensional space, the spin breaking terms can only be in the form of a ■ cr, where a is an arbitrary three 

dimensional vector. After analyzing all the possible combinations of the field and & operators, we find that there are 
three independent spin-breaking terms at sub-leading order in 1/itiq, which are given by 

-4s = Sfs^iW^fo-^E 1 + h.c. , (10a) 

£ S Q P L = -iJ^—ejikfjrijio-jjk] - Tr[7V7- , '])£'' + h.c. , (10b) 
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Table 1 : The numerical values of the coupling constants S j (GeV 1 ) determined by fitting the experimental decay widths. 



Charmonium 


Bottomonium 


Decay Width (keV) 


^'^(GeV- 1 ) 


Decay Width (keV) 




rOfco -» J /iff + y)= 121.7 ±10.9 


(2.13 ± 0.95) X 10" 1 


r(T(2S)^Xbo(lP) + r) = 1-22 ±0.16 


(9.01 ± 0.50) x 10" 7 


rCfci -» //^ + y) = 295.8 ± 21.5 


(2.31 ± 0.08) X 10" 1 


r(T(25) -> XbiO-P) + T) = 2.21 ± 0.22 


(9.89 ± 0.50) x lO" 7 


T(Xc2 — * //^ + y) = 384.2 + 26.6 


(2.29 ± 0.08) x 10-' 


r(T(25) -» Xb2(lP) + y) = 2.29 ± 0.22 


(9.86 ± 0.47) x 10-* 


r(f->AT c o + y) = 29.2 ±1.3 


(2.25 ± 0.04) x 10"' 


r(T(3S)^Xbo(2P)+y) = 1.20 ±0.16 


(1.39 ± 0.09) x 10" 1 


iW-># c i + y) = 28.o ±1.5 


(2.36 ± 0.06) x 10-' 


T(Y(3S ) -» ^ M (2P) + r) = 2.56 ± 0.34 


(1.57 ± 0.10) x 10" 1 


W -»AT«fl+y) = 26.6 ±1.3 


(2.74 ± 0.07) x 10"' 


r(T(35) -» XbidP) +y) = 2.66 ± 0.41 


(1.55 ± 0.12) X 10- 1 



£ SP T = -X_(Tr[/V^V] + Tr[to-7V 7 ])£' + h.c. , (10c) 



znP,mS 



-QT 



where S'fi , S n /f lS and S'l p ' mS are the coupling constants that are suppressed by 1 /m?,, since the spin-breaking poten- 



tials is of 0(1 /w^) compared to the static potential! 15 1 



After including the spin-breaking contribution, the formula of the El transition decay widths turn to be: 

t i {°oj > % M„ P 

T(m 3 S l _ n 3p / ) = (2J + l)-Sy lla) 

9tt y M mS 

/ cfiPjnS \2 

r(n -> m Si) = — ^ y T^— , (Hb) 

37T 7 M nP 

T^S^n^,) = ^—A (11c) 

/ cnP,mS \2 

r(« Pj -> to S ) = — k —— , (lid) 

3n 7 M„ P 

where k y is the energy of the emitted photon, and 

cnP,mS cnP,mS cnP,mS , o cnP,mS , o cnP,mS / 1 o \ 

°Q0 = 6 Q ~ 5 QS +26 QL +3 V (12a) 

cnPjnS ciiPjnS ctiP,mS , ctiP.mS ^ ctiP.mS /1 
6 Gl = (5 Q _<5 GS +5 Gi _2 V ( 12b ) 

S Qi mS = ^ ~ 6 QS mS - 6 "qL S + X 6 QT mS ( 1 2C > 

ctiP.mS s-nP.mS , o c-nP.mS n w etiP.mS , c-nP.mS /i <-»j\ 

In principle, t^y'" 5 can be obtained by calculating the matrix element of the electromagnetic current between the 
wave functions of the n 3 P j and m 3 S i states in pNRQCD (or in any potential model, see Ref. |0] for a recent review). 
For some processes, such as XcJ — * JNf + y> - > XcJ + 7» Y(2S) — > ^w(lP) + y, and T(3S) — > Xbj(2P) + y, their 
decay widths have been measured [16], so we can obtain the values of the corresponding coupling constants S'gj nS 
(for J=0, 1 ,2) by fitting the data, which are list in Table 1 . Note, in our treatment of the uncertainties, we only take into 
account the uncertainties in total decay widths and the branching ratios. The formulas in Eq. dT2b show that up to the 
sub-leading order tfgf (for J=0,l,2) only depends on 5q MS - 6"^, 6"^ s , and 8"^. After resolving Eq.(Q2}, we 
obtain that 

glP.VI _ 6 IPIS = 22 J x 1Q -2 GeV -l 6 lP,2S _ 6 IP2S = 6 x l0 -2 Ge y-l f 



J cS " wtv '"c "cS 

S\ p ' ls - 5 l b f s = 9.7 x lO- 2 GeV-\8f- 3S - d^' 25 = 15.3 x lO^GeV' 1 , (13) 



Table 2: The values of the coupling constants for the spin-breaking terms S"q^" and $qt determined by fitting the experimental decay widths 
(unitlO^GeV- 1 ). 



Ap.lS 
°cL 


o cT 


a p ;is 

°cL 


"cT 


Ma 

°hT. 


a p ;is 

°hT 


°hT. 




-0.3 ± 0.4 


-0.4 + 0.3 


-1.8 ±0.3 


0.1 ±0.2 


-0.1 ±0.3 


-0.1 ±0.2 


-0.2 ± 0.6 


-0.3 ± 0.3 



The values of the corresponding 6"™, and 5 nP ^ lS are given in Table 2. If we assume that 6 nP f s , 6 nP [ nS , and 5 nP ^ ,s 
are in the same order, the results in Table 2 and those in Eq.([r3l will indicate that in these decay processes the 
contribution of the spin-breaking effect is less than that of the leading order term by at least a factor of 10, and 
furthermore comparing to the 25 — > IP transition process in charmonium system the spin-breaking effect in the 
bottomonium system 25 — » IP process is of (m c /mh) 2 suppressed, which is consistent with the power counting rule 
of pNROCD llll. 



3. Radiative Transitions of if/(l 3 D j>) -> Xcj + 7 

The spectrum of Z)-wave heavy quarkonia has been calculated in potential model by many groups, fox example 
recently in Ref . Jl7l [l8ll . In cc system, the (^(3770) state is treated as a pure 3 D\ state or a predominant D-wave state 
with a small admixture of 25 state |T9l 2(J. The other states 3 £>2, 3 £>3 and X D^, whose decay widths are all expected 
to be narrow, have not been observed yet. The masses of tfr( 3 D2) and ria^D?) predicted by potential models lie 
between DD and DD* thresholds [17]. They are forbidden to decay into pseudoscalar pair by parity. The narrowness 
of tfr^Dj,) is due to that its decay into DD is a P-wave decay, which is highly suppressed. Hence, the branching 
ratios of their radiative decays are all considerable. The widths of the radiative transition nD to mP can be calculated 
straightforwardly by employing the Lagrangian in Eq.dTbl: 



/ cnDjnP\2 



in 



M, 



(14a) 



nD 



r(n'D 2 -> m'Pi) = 



/ cnD,mP\2 

( 6 Q ) l3 M mP 



(14b) 



where the coefficients S j>j are S ij = 5/9, 5/12, 1 /36, S2,j = 0, 3/4, 1 /4, and S 3 y = 0, 0, 1 for J — 0, Irrespectively 



The decay widths of i^(3770) decay toxcQ + J and^ c i + y given in PDG are lfl6ll : 

T(^(3770) -> Xcoj) = 199 ± 26 keV, T(i/r(3770) ^ X c\j) = 79 ± 17keV. 



(15) 



If (^(3770) is a pure ID state, the values of the coupling constant 6 l c Dlp determined through T(i^(3770) — > Xcj + J 
are 6 l c Dlp = 0.31 ± 0.02 GeV" 1 , and 5\ mp = 0.35 ± 0.04 GeV" 1 for J = and J = 1, respectively, which are very 
close to each other. In this case, it indicates that the heavy quark spin symmetry is roughly preserved in the ID to IP 
radiative transitions. The average value is 



~5 l c Dlp = 0.32 ±0.02 GeV" 1 . 
If we treat (/r(3770) as a S-D mixing state and using the same notation in Ref. I20I1 : 

i/<3686) = cos(6»)[25> - sin(0)|l£>>, (A(3770) = cos(0)|l£>> + sin(6»)|25), 



(16) 



(17) 



1 For the 3 I>i decay into 3 Pj + y, our results do not agree with those in Ref.| 14]- After private communications, their new results in the erratum 
(2lll now agree with ours. 
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the analytical formulas for i^(3686) and i/f(3770) decay into y + Xcj turn to be 



1^(3686) -> X c0 +J)= a J" k l (cos 2 (0)(c5^ 2s ) 2 - 2 J- sm(0) cos(0)^ p - 2 ^' D - lp + - sin 2 (0)(^ D ' 1P ) 2 ) (18a) 



3-^r e i jJS , 2rn\rxlP2S\2 , ^ ^;„/-a\ „ n cra\xWS X 1D,1P , $ c,;„2/ Q ^ADAP^ 



r(^(3686) ^^ cl +r) = *" ^(cos z (fl)(tf'^) z + J - sin(0)cos(0)(5r'"<5r ir + sin^X^'T) (18b) 

y7ryw^ (36 86) y 3 12 

r(^(3686) -^Xc2 + 7) = 5Mfe ^(cos 2 ^)^' 25 ) 2 - J-^ sin(6»)cos(0)^ 2S (J ( 1A1/> + J- sin 2 (0)(^ D ' 1F ) 2 ) (18c) 



V(3686) 



15 v ' v ' c e 60 



r(<A(3770) -» * c o + r) = -j 5 ,?* 4 (cos 2 (0)(^ D ' 1F ) 2 + 2 J| sin(0) cos^)^' 1 ^ 25 + | sin 2 (0)(<^ 2s ) 2 ) ( 1 8d ) 



5 v 7 v y e c 5 



5M /T 12 

r(^(3770) = . A ./" (cos 2 (g)(^ lf ) 2 -4 J - sin(g) cos(g)^ D - lp ^- 2s + — sin 2 (fl)(c^ 2s ) 2 ) (18e) 

367TAfy(3770) V 5 5 

T(^(3770) ^^ c2 +r) = ino A ff ^ (cos 2 (0)(^ D ' IP ) 2 +4 VTSsin^cos^^'^^+eOsin 2 ^)^ 25 ) 2 ) (18f) 

The mixing angle = (12 + 2)° that is determined from the leptonic decays of i/f(3686) and i/f(3770) also is favored 
by some other considerations [22]. If we fit i/f(3770) and (i/f(3686)) decays into^o + J with = 12°, we obtain two 
set solutions, which are labeled by subscripts 1 and 2, respectively, 

d\ p ' 2s = 0.30GeV~ I ,<5| D ' lp = 0.27GeV _1 ; S l 2 F ' 2S = 0.14GeV" I ,^ D ' IP = -0.34GeV _1 ; (19) 

If we fit their decay into^ c j + y the results are: 

S\ P ' 2S = 0.18GeV _I ,£j EUP = 0.42GeV _1 ; 6 l 2 P ' 2S = 0.28GeV" I ,^ D ' IP = -0.27GeV -1 ; (20) 

The difference between the results in Eq.([T9l> and those in Eq.(l20b shows that the heavy quark spin symmetry is 
largely violated if i/f(3686) and (^(3770) are assigned as two S-D mixing sates as given in Eq.(fT7b with mixing angel 
9 — (12 + 2)°. Consequently, to understand the radiative decays of </f(3770) and </f(3686) in the S-D mixing picture, 
some other effects like the relativistic corrections[23] or the couple channels effect[24] should be taken into account. 

Since in the S-D mixing picture the heavy quark spin symmetry does not hold anymore, thereafter we will adopt 
that i/f(3770) is a pure D-wave state and choose the value of the coupling constant to be that in Eq.([T6ii to study the 
radiative decay of the D-wave states. The upper limit of i/f(3770) — > Xc2 + y is that S(t/'(3770) — » \c2 + y) < 9 x 10~ 4 
Ills]. Using the result in Eq. (ll4ab . we predict that 

T(i/r(3770) -> Xc2 + 7) = 2.55 ± 0.28 keV, £(0<377O) -> Xo2 + y) = (9.4 ± 1.0) x 10~ 5 . (21) 

which is compatible with the experimental data and is about 4 times larger than those in Ref.[ 14]. The i/r(3770) decay 
into Ycj + 7 nas a ls° been studied by potential model. For comparison, we choose two potential models calculations 
lil9[ll7ll . in which the predictions of i/f(3770) decay into^f f o,i + y agree well with the experimental data after including 
relativistic corrections. Their predictions of r(i/f(3770) — > Xdi + 7) ws 3 .0 01911 or 3 .3 111711 keV. Both of them are 
consistent with our results. 

As mentioned above, the other ID states are all expected to be narrow. Their spectrum and the El transition decay 
widths have also been calculated in Ref ill 711 . Their results are: 

M(1 3 D 2 ) = 3.838GeV, T(1 3 D 2 -> Xd(Xc2) +y)= 268(66) keV (22a) 

M(1 3 D 3 ) = 3.849GeV, T(1 3 D 3 -> Xc2 +j) = 296 keV (22b) 
M(1 1 D 2 ) = 3.837GeV, T(l 'D 2 -> h c + y) = 344 keV. (22c) 
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If we choose the same mass values, our predictions are 



r(l 3 D 2 -> Xc\ +j) = (288 ± 25) keV, T(1 3 D 2 -> xa +j) = (50.3 ± 5.5) keV, (23a) 

r(l 3 D 3 -> Xci + r) = 224 ± 25 keV, r(l 1 D 2 -> h c + y) = 267 + 29 keV. (23b) 

which agree with the potential model results. 

Recently, the X(3872) state has received much attention since it was first discovered by Belle Collaboration 112511 . 
and then was confirmed in pp collision at Tevatron 1 26] . It was also observed by Babar Collaboration Until now, 
there is not a convincing explanation about its nature yet. Only the charge parity C = + is established from its decay 
into J/ip + y JUt]. After analyzing B — > J/ip + oj + K, Babar Collaboration found its J PC favors 2~ + lE9ll . If it is a pure 
charmonium D-wave state, the only assignment will be the J] C 2( l D2). It then should has a sizeable decay into y + h c . 
Evaluating in a similar way, we obtain 

T(X(3872) -> h c + y) = 359 + 59keV, (24) 

which is very large. So studying X(3872) decay into y + h c will be helpful to understand its nature. 

The D-wave bottomonium states were observed from the cascade of T(35) [30], however no other further in- 
formation is known yet. We can not make any prediction about their radiative decay with the effective Lagrangian 
method at present. 



4. Relation Between Radiative Decay and LH Decay of Xbj(nP) 

The total decay widths of the P-wave bottomonium states XbAnP) (n = 1,2) have not been measured yet, so 
we can not compute S'f' by fitting the data. Besides the radiative decay, the LH decay is also an important decay 
mode for the P-wave quarkonium states. One remarkable success of NRQCD is that it can systematically resolve the 
infrared divergence problem in the CS model (CSM) calculation for the LH decays of P-wave states by introducing 
the CO contribution H, 31]. For the states in strong coupling region, where most of the heavy quarkonium states 



below threshold are expected to belong to, further study of pNRQCD shows that the CO matrix elements can be 
related to the wave function of the bound states llOlll ill . In particular, in the strong coupling region the ratio pAnP) = 
mfH&(nP) fHiinP) does not dependent very much on the radial quantum number n 11 Oil . By fitting the open charm 



decays of>w(lP) and*w(2P), CLEO collaboration obtained thatp g (lP) = 0.160+°™ andp 8 (2P) = 0.074+°°'° J32J 



which is a little different from pNRQCD prediction. 

Next, we will show that the relation between the radiative decay and the LH decay could provide another way to 
extract the values of p%{\, 2P). According to NRQCD approach, at v/, leading order, the LH decay width for P-wave 
states is given by: 

r(Xbj(nP) -> LH) = — + — (25) 

m b m b 

where C \ and Cg are the /-dependent short distance coefficients and have been calculated up to a- 3 order [ 34] . Although 
neither of the radiative and LH decay widths have been measured, their branching ratios are known. Recently Babar 
Collaboration update the branching ratios of Xbj(nP) — > T(m5) + y, their latest results are j33ll : 

SxbAlP) -> T(1S) + y = (2.2 ± 1.5^;° ± 0.2, 34.9 ± 0.8 ± 2.2 ± 2.0, 19.5 ± 0.7!j; 3 ± 1.0)% for (J = 0, 1,2) (26a) 
&Xbj(2P) -> T(25) + y = (-4.7 + 2.8^ ± 0.5, 18.9 ± 1.1 ± 1.2 ± 1.8, 8.3 ± 0.8 ± 0.6 ± 1.0)% for (J = 0, 1,2) (26b) 
&XbA2P) T(15) + y = (0.7 ± 0.4^ ± 0.1,9.9 ± 0.3^4 ± 0.9,7.0 + 0.2 ± 0.3 ± 0.9)% for (J = 0, 1,2) (26c) 

The branching ratio of XbAnP) decay into LH can be obtained by subtracting its all the known transitions to other 
bottomonium states, which can be read out directly from PDG llol . The ratio of the two branching ratios can be 
expressed as 

j? , m ^XbAnP) -» LH) T(x b AnP) -» LH) 3nM„ P n + rrvut ™ 

RAnP) = = = -(CAuYnAnP) + C%{uyrl%(nP)) 

S(x b AnP)^y + r(nS)) r(x b AnP)^y + T(nS)) M„ s 6" p '" s k] V lW 8W 

(27) 
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Rj(nP) only dependents on three unknown parameters "HxinP), "H%{nP) and 6" b P '" s . Therefore, to compute the ratio 
p%{nP), we only need two independent inputs. Since the uncertainties of th&Xbo(l,2P) radiative decays are are large, 
we choose the data of Xb\,i(^,2P) decays. Using the a 3 s order short-distances coefficients listed in Ref. 13411 and 
setting — fi — 2mt, a s (2mi,) = 0.18, and the number of light flavor quark Nj = 4, we obtain p&(lP) = 0.150^qq^ 
and ps(2P) = 0.110 + 0.030. Our value of pa(lP) is a little smaller than that of CLEO, while our value of p s (2P) 
is about 1.5 times larger than that of CLEO, which makes p$(lP) close to p%(2P). This indicates that the pNRQCD 
assumptions is reasonable to study the LH decays of xui^P) and^/,y(2P) states 11 Oil . 



5. Summary and Conclusion 

In summary, the radiative decays of the heavy quarkonia are studied with the help of the effective Lagrangian. 
To have a better understanding of the radiative transitions among S - and P— wave states, we take into account the 
contribution that is due to the spin-breaking interactions. By fitting the experimental data, the coupling constants of 
the spin-breaking terms in i//{2S), XcJ< Y(2S) and T(35) radiative are obtained, which are listed in Table 2. We find 
that the values of the coupling constants in the spin-breaking terms are less than those in the leading order terms by at 
least a factor of 10 and that the spin-breaking terms in i//(2S ) and T(2S ) indicate that the spin-breaking contribution is 
suppressed by l/m^, which agrees with pNRQCD power counting rule. We also calculate the radiative decays of the 
cc(lD) states, whose total decay widths are expected to be narrow. Based on that i/r(3770) is a pure D-wave state, our 
predictions of the radiative decay widths of the other D-wave states are consistent with the potential model results. 
Furthermore, we predict that T(X(3872) -> h c + y) = 359 + 59keV, if X(3872) is the 2~ + state. We also study the 
S-D mixing effect in 1^(3770) radiative decay and find that there is no heavy quark spin-symmetry if the mixing angel 
is 12°. As an useful application, we find that relating the radiative decay of Xbj(nP) to their LH decays can provide 
another way to estimate the ratios of p%{nP) - m 2 h ( H%(nP)l < H\(nP). By fitting the data of Xbii^P, 2P) and that of 
Xb2(lP, 2P), we get that p 8 (lP) = 0.150+^ and Ps(2P) = 0.1 10 + 0.030, which approximately equal to each other. 
Our result provide an evidence on pNRQCD assumptions lioll . 
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